Resolution in Focused Electron- and Ion-Beam Induced Chemical Vapor Deposition 
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The key physical processes governing resolution of focused electron-beam and ion-beam-assisted 
chemical vapor deposition are analyzed via an adsorption rate model. We quantify for the first time 
how the balance of molecule depletion and replenishment determines the resolution inside the locally 
irradiated area. Scaling laws are derived relating the resolution of the deposits to molecule dissocia- 
tion, surface diffusion, adsorption, and desorption. Supporting results from deposition experiments 
with a copper metalorganic precursor gas on a silicon substrate are presented and discussed. 
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Local electron-beam induced chemical vapor deposi- 
tion (CVD) is a physical phenomenon well known from 
the build up of carbon contamination since the begin- 
ning of electron microscopy Recent research has 
spent tremendous efforts on the systematic creation of 
functional nanoscale deposits by means of focused elec- 
tron beams and - with the development of scanning ion 
microscopes - of focused ion beams. Organic, organo- 
metallic, and inorganic precursor gas molecules were sup- 
plied into the microscope chamber. Upon irradiation, 
deposition results from non-volatile dissociation prod- 
ucts whereas etching occurs when a reaction of disso- 
ciation products with the substrate leads to the forma- 
tion of volatile species. These local deposition and etch- 
ing techniques have numerous potential applications in 
nanosciences including fabrication of attachments in me- 
chanics , high-resolution sensors in magnetic, thermal, 
and optical scanning probe microscopy d, 0, Q , optical 
elements in nanooptics @, 0] , contacts in electronics Q , 
and nanopores for ionic current measurements of cells 
and DNA in biology @,[I3|- 

Few experiments analyze the physics of focused 
electron-beam (FEB) and focused ion-beam (FIB) in- 
duced deposition and etching. For FEB, Allen et al. 
pd| noted that the spatial flux distribution of electrons 
passing through the adsorbed molecule layer consists of 
the incident primary beam (of up to several keV) and 
the emitted electrons: back scattered primaries and low- 
energy (< 50 eV) secondary electrons. The entire energy 
spectrum is responsible for the dissociation process. For 
FIB, Dubner et al. [ijj showed that this spectrum is asso- 
ciated to the energy deposited into the substrate surface 
through the collision cascade generated by the primary 
ion beam. 

Our work is motivated by the fact that focused parti- 
cle beam induced chemical vapor deposition and etching 
processes are widely used in nanoscale fabrication, but 
there are only very few attempts to describe the spatial 
resolution of this process theoretically. For a singular pri- 



mary electron beam Silvis-Cividjian et al. [H| and Ha- 
gen et al. [T3| concluded from Monte Carlo simulations 
that the ultimate resolution depends on the emitted sec- 
ondary electron distribution. However, their assumption 
that the irradiated area is permanently covered with a 
monolayer of adsorbed molecules is idealized. In fact, 
as we will show, this coverage results from a balance of 
molecule depletion by dissociation and molecule replen- 
ishment strongly depending on adsorption, desorption, 
and diffusion. This is also White's et al. conclusion [l5| 
who studied the gas transport phenomena at the micro- 
scope chamber scale affecting the overall deposition rate. 
Midler's model 16] for FEB induced deposition, which 
was later adapted by Haraichi et al. [12] to gas-assisted 
ion-beam induced etching, takes the key processes of sur- 
face diffusion, dissociation, desorption, and adsorption 
into account via an adsorption rate equation. However, 
he used a flat-top beam distribution which allows no con- 
clusions about resolution. 

In this letter we present an adsorption rate model with- 
out the above mentioned limitations, considering two rel- 
evant peak distributions for the incident beam and for 
the emitted spectrum, and a non-dissociative Langmuir 
adsorption term. It allows to derive two scaling laws for 
resolution and to estimate important physical parame- 
ters of the process. Finally, we discuss results of carefully 
designed experiments, clearly supporting the theoretical 
conclusions. 

The model assumes second-order kinetics of molecule 
dissociation by electrons. In a system with rotational 
symmetry the vertical FEB deposition or etch rate R{r) 
(in units of dimension per unit time) as function of the 
distance r from the centre of the primary electron (PE) 
beam is thus: 



R(r) = Vn{r) 



a(E)f(E,r)dE^Vn(r)af(r) (1) 
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where V is the volume of the decomposed molecule or 
etched atom, n(r) is the number of adsorbed molecules 
per surface unit, a{E) is the energy dependent electron 
impact dissociation cross section, Epe is the energy of 
the PEs, and f(E, r) describes the spatial flux distri- 
bution of the electron energy spectrum generated by 
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FIG. 1: Reference system and schematics of processes in- 
volved in FEB induced deposition. Inside the irradiated area 
precursor molecules are depleted by dissociation and replen- 
ished by adsorption and surface diffusion (dashed arrows). 
Symbols are defined in text. 

the PEs. Since the energy integral can be solved only 
approximatively due to missing a(E) data of adsorbed 
molecules and uncertain parameter estimates entering in 
the Monte Carlo simulation of the emitted electron dis- 
tribution [Til ] . we use the simplified expression in Eq. 
[TJ where a represents an integrated value over the en- 
ergy spectrum. Such cross sections were measured for 
several relevant molecules [H [U HH, [H, [H, The 
spatial distribution /(r) is a convolution of the Gaus- 
sian incident beam distribution f(r) oc exp(~r 2 ) with 
a peak function for the emitted spectrum which can be 
approximated by f(r) oc exp(-r). Full widths at half 
maximum (FWHM) of the emitted distributions range 
between - 0.1 nm (200 keV) and 2 nm (1 keV) 
Similar considerations apply for FIB where the spatial 
distribution being responsible for molecule dissociation 
is a convolution of the primary beam distribution with 
the distribution of excited surface atoms generated by 
the collision cascade [TtJ • Table UJ summarizes typical 
FWHM values of primary Gaussian beams. 

Four key processes as depicted in Fig. Q] are consid- 
ered to determine the surface density n(r) of adsorbed 
molecules: a) adsorption from the gas phase governed 
by the precursor flux J, the sticking probability s, and 
coverage n/no] b) surface diffusion from the surround- 
ing area to the irradiated area governed by the diffusion 
coefficient D and the concentration gradient; c) desorp- 
tion of physisorbed molecules after a residence time t; 
d) molecule dissociation governed by the product af{r). 
For the molecule adsorption rate dn/ dt follows: 

dn j ( i n \ £> ( ® 2fl ^ dn\ 
dt \ n J \ dr 2 r dr J 

Adsorption Diffusion 
77 

- afn .(2) 

Decomposition 

Dcsorption 

The adsorption term in Eq. [^describes a non-dissociative 



Langmuir adsorption, where no is the maximum mono- 
layer density given by the inverse of the molecule size. 
This adsorption type accounts for surface sites already 
occupied by non-dissociated precursor molecules and lim- 
its the coverage to n$. All parameters other than n — 
n(r,t) and / = f(r,t) are considered constant. 

Solving eq. [2] for steady-state (dn/dt = 0) and neglect- 
ing the diffusion term we obtain n(r) — sJT c ff(r) with 
the effective residence time of the molecules T c s(r) — 
(sJ/no + 1/r + <rf(r)) . The deposition or etch rate 
becomes 



R(r) = sJT eS (r)Vaf(r) (3) 

and represents the deposit or etch shape at a given 
time. For any peak function f(r) with a peak value 
fo = f(r = 0), we can define the effective residence 
time in the center of the eletron beam r; n = r o fr(0) = 
l/(sJ/n + 1/t + afo) and the effective residence time 
far away from the electron beam center r ou t = T c g(r — > 
oc) = l/(sJ/rto + The dimensionless ratio f = 

^out/Tin = 1 + o"/o/(l/r + sJ/no) represents a mea- 
sure for depletion of precursor molecules due to dissoci- 
ation at the center of the beam. With a Gaussian beam 
f(r) = fo exp(— r 2 /2a 2 ) we derive the first scaling law of 
deposit size as function of depletion: 

<p = {log 2 (l + f)} 1/2 = FWHM d /FWHMb, (4) 

where FWHM B and FWHM D are the full widths at 
half maximum of f(r) and R(r) or, in other words, the 
FWHMs of the incident beam and the deposit. The 
idealized case of zero depletion, i.e. permanent mono- 
layer coverage, corresponds to f = 1. Then deposition 
or etching proceeds in the electron-limited regime and 
the deposit (or etch) size corresponds to the electron 
beam distribution since the logarithmic term becomes 
1. With increasing depletion the deposit (or etch) size 
becomes steadily larger than the beam size. For the 
peak function /(r) oc exp(-r) the expononent in eq. 0] 
becomes 1. The degree of depletion strongly depends 
on the dissociation frequencies afo summarized for FEB 
and FIB in Table [TJ In order to replenish the dissoci- 
ated molecules inside the continuously irradiated area by 
gas transport only, we need f — > 1, i.e. the precursor 
supply frequency sJ/no should exceed cr/o, being equiv- 
alent to > 2 x 10 3 ML/s (monolayers per second) for FEB. 
This corresponds to a precursor flux on the substrate of 
J = 2xl0 17 molecules cm~ 2 s _1 , settings = 1 and taking 
no = 10 14 cm~ 2 as typical value. For FIB several orders 
of magnitude larger gas supply would be needed. Desorp- 
tion frequencies are situated around t -1 = 10° . . . 10 3 Hz 
[Tel , [l9| . Above estimations clearly suggest that most of 
the FEB and FIB processing experiments were performed 
in the precursor-limited regime limiting the minimum de- 
posit or etch size. 
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TABLE I: Typical ranges of incident peak flux /o and size 
FWHMb of a focused electron beam (5 keV, field emission) 
and an ion beam (30kV, Ga + ). Representative ranees for a 
were collected from [H, EH |Hj for FEB and from HI [H, HI 
for FIB. 





fo 


FWHMb 


a 






[l/nm 2 s] 


H 


[nm 2 ] 


m 


FEB 


8 x 10 6 . . . 


2.5. . . 


2 x 10" 4 . . . 


2 x 10 3 . . . 




5 x 10 7 


100 


2 x 10" 1 


1 x 10 7 


FIB 


2 x 10 5 . . . 


7 ... 


10. . . 


2 x 10 6 . . . 




5 x 10 6 


100 


50 


2.5 x 10 8 



In the following we show under which conditions sub- 
stantial replenishment by surface diffusion can be ex- 
pected. Solving equation [5] numerically with MATLAB 
for steady-state (dn/dt = 0), the boundary conditions 
n(r — * oo) = n out = sJr out and dn(r — 0)/dr = 0, and 
with a Gaussian distribution f(r) = f cxp(— r 2 /2a 2 ) we 
get n(r) and finally R(r). A plot against the dimension- 
less variable f — 2r/FWHMB for a given depletion and 
diffusive replenishment results in a unified representation 
of deposit shapes for any FWHMb , see Fig. [2j Replen- 
ishment by diffusion is described by the dimensionless 
ratio p = 2pi n /FWHMB, relating the diffusion path in 
the center of the beam p m — {Dt^) 1 ! 2 with D being 
the diffusion coefficient, to the beam size. For p = 0, 
the flat top shape is defined by Eq. [3] With increas- 
ing diffusive replenishment deposits change into indented 
and round apex shapes, since adsorbed molecules increas- 
ingly reach the centre of the irradiated area before being 
dissociated. Hence both deposition rate and resolution 
increase. The maximum diffusion enhancement in depo- 
sition rate becomes R(p = oo) / R(p = 0) = f at r = 0. 
For p — > oo, Eq. [3] simplifies to R(r) = n out Vaf(r) since 
any depletion is entirely compensated by diffusion and a 
permanent monolayer coverage provided. In other words, 
the electron-limited regime is established and the deposit 
shape corresponds to the electron beam distribution /(r). 

Figure [3] represents a graph relating the dimensione- 
less deposit size (p = FWHMd/FWHMb to irradiative 
depletion and diffusive replenishment. For any depletion 
with p > 2 we get Cp < 1.03, i.e. the deposit size is within 
3% close to FWHMb when the related diffusion path in- 
side the irradiated area becomes at least comparable to 
the size of the electron beam distribution. The deposits 
become broader than the electron beam for f > 1 and 
small p, branching out into constant maximum size given 
by Eq. U at negligible diffusive replenishment p — > 0. 
Figure 3 holds independently of how diffusive replenish- 
ment is experimentally achieved: either via the beam size 
FWHMb (using the focus of the beam) or via the diffu- 
sion path p ln (changing precursor diffusion) . The second 
scaling law of deposit size as function of diffusive replen- 
ishment is obtained as (circles in figure [3]) : 
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FIG. 2: Normalised steady-state deposition rate at indicated 
depletion from Eq. [1] representing the deposit shape. The 
diffusive replenishment p = 2pi n /FWHMb is varied. Note 
the shape transition from flat top, p — 0, indented, p = 0.32, 
to Gaussian, p = oo. 
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FIG. 3: Normalized deposit size vs. normalized diffusion path 
for varying depletion (indicated). At p — 2 the diffusion path 
equals the beam size. The curve for f = 10.9 corresponds to 
the FWHMs of shapes in Fig. [2] Circles represent Eq. [5] The 
inset shows the FWHMd definition of indented deposits. 



0«{log 2 (2 + /T 2 )} 



1/2 



(5) 



For the exponential peak function f(r) oc exp(-r) the 
expononent in Eq. [5] becomes again 1. Both Eq. [4] and 
Eq. Ogive upper limits for <p(f) and <p(p). The smaller 
value of both defines the minimum deposit (or etch) size 
with respect to the beam size. 

Now we present FEB deposit shape measurements 
obtained with a special designed SEM compatible 
atomic force microscope (AFM). The advantage is that 
AFM reveals three-dimensional topography features, 
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FIG. 4: AFM image and line scans of FEB deposits from 
Cu(hfac)2 precursor. Exposure times are indicated. The in- 
dented apex shapes are due to depletion. Dashed lines repre- 
sent fits of Eqn. [T] 

like indented shapes, which are difficult to resolve in 
SEMs due to edge contrast effects. We used Cu(II)- 
hexafluoroacetylacetonate precursor molecules impinging 
on a native Silicon substrate with J /no = lOML/s un- 
der irradiation with a 5 keV Gaussian electron beam 
(jo = 9 x 10 4 iim-V 1 and FWHM B = 110 nm). An 
indented shape with FWHMd = 200 nm is observed, 
hence (p = 1.8. Assuming r = 10~ 3 s and s — 1 we 
get from Eq. 2] f = 8.9, i.e. a ~ 0.09 nm 2 . From 
Eq. [5] follows p — 0.37, i.e. p ln — 20 nm. Using 
the relation p m ~ (D/afo) 1 ^ 2 for f 3> 1 results in 
fl~3x 10~ 8 cm 2 s -1 . 

The values for the depletion, cross section, and diffu- 
sion coefficient represent lower limit estimations since the 
same FWHM-ratio (p can be obtained with larger deple- 
tion and larger diffusive replenishment, see figure[3] Tak- 



ing as maximum dissociation cross section the molecule 
size a = 0.6 nm 2 [25j |. we get f = 60. The corresponding 
upper limit estimate for the diffusion coefficient is derived 
from the shape fit in Fig. 2] to be D = 4 x 10 _7 cm 2 s _1 . 
Finally, we estimate diffusion coefficients needed for es- 
tablishing the electron-limited regime under typical irra- 
diation conditions summarized in table HI Compensation 
of depletion by surface diffusion requires a beam focus 
smaller than the molecule diffusion path, p > 2, see fig- 
ure El Together with the relation p- m ~ (D/afo) 1 / 2 for 
f » 1 we obtain D = 10" 12 . . . 10~ 6 cm 2 s" 1 for FEB 
and D = 10~ 7 . . . 10~ 3 cm 2 s" 1 for FIB. 

In conclusion, we quantified the crucial role of deple- 
tion and diffusive replenishment on deposit resolution for 
two relevant distributions: an incident Gaussian beam 
and an emitted distribution with exponential decay via 
an adsorption rate model. Our model is applicable to 
gas-assisted deposition and etching with focused electron- 
and ion beams and thus covers numerous applications 
in fundamental research and in nano-scale fabrication. 
We demonstrated how physical parameters can be de- 
termined from fitting experimental deposit shapes with 
our model. An extension of the studies to different beam 
shapes or experimental arrangements like pulsed beams 
at different temperatures is straightforward. It will en- 
able the systematic determination of all physical key pa- 
rameters involved in the process thus opening the door 
to the controlled fabrication of tailored nanoscale devices 
by charged particle beam induced CVD and etching. 

We acknowledge financial support from the European 
Commission, FP6 Integrated Project NanoHand (IST-5- 
034274). 
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